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ABSTRACT 

Observational evidence for black hole spin down has been found in the nor- 
malized light curves of long GRBs in the BATSE catalogue. Over the duration 
Tg of the burst, matter swept up by the central black hole is susceptible to non- 
axisymmetries producing gravitational radiation with a negative chirp. A time 
sliced matched filtering method is introduced to capture phase-coherence on in- 
termediate timescales, r, here tested by injection of templates into experimental 
strain noise, h n (t). For TAMA 300, h n (f) ~ 10" 21 Hz~t at / = 1 kHz gives 
a sensitivity distance for a reasonably accurate extraction of the trajectory in 
the time frequency domain of about D ~ 0.07 — 0.10 Mpc for spin fown of black 
holes of mass M = 10 — 12M Q with r = 1 s. Extrapolation to advanced detectors 
implies D ~ 35 — 50 Mpc for h n (f) ~ 2 x 10~ 24 Hz~2 around 1 kHz, which will 
open a new window to rigorous calorimetry on Kerr black holes. 



1. Introduction 



Calorimetry and spectro scopy on all radiation channe ls will be key to identifying Kerr 
black holes in the Universe ( Ivan Putten fc Levinsonl l2002l ). They may be central to some 
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of the core-collapse supernovae (CC-SNe; IWoosleyl (119931 )) and mergers of ne utron stars 



with a nother neutron star (NS-NS) or with a companion black hole hole (BH-NS; iPaczyhski 



( Il99ll ))-the currently favored astronomical progenitors of cosmo logical gamma-ray bursts 
(GRBs). While all short GRB are probably produced by mergers, the converse need not 
hold in general. B lack holes with a neutron star companion should have a diversity in spin 
( Ivan Puttenl Il999l ). These mixed binaries can be short and long lived, depending on the an- 
gular velocity of the black hole. It predicts that some of the short GRBs , produced by slowly 



spinning black holes , also feature X-ra y afterglows (Ivan Putten 



(IGehrels et al.l 120051 ) and GRB050709 (IVillasenor et all 120051 iFox et all 120051 iHiorth et al 



200 LaJ) as in GRB050509B 



20051 ). Long GRBs, produced by rapidly r otating black hole s, are expected to fo rm in CC- 
SNe from short, intra-day period binaries (iPaczvriskil (ll998f); Ivan Puttenl (120041 )). but also 
from mergers with a companion neutron star ( van Puttenl Il999l ) or out of the merger of two 



neutron stars (IBaiotti et all 120081 : Ivan Putten 



van Putten 
20091). 



A diversity in the origin of long GRBs in CC-SNe and mergers naturally accounts for 



event s with and without supernovae, notably GRB0606 14 (Ivan Putten 



2008 



2009 ). with and without pronounced X-ray afterglows ( van Putten fc Gupta 20091 ) and in 



Caito et al 



wind versus constant density host enviro nments, that are re levant to recent studies of ex- 
traordinary Swift and Fermi-L AT events ( Cenko et al 2010al lb[). 



Rapidly rotating black holes can sweep up surrounding matter and induce the formation 
of multipole mass-moments. In this process, spin energy is cat a lytically converted to a 
l ong d uration gravitational wave burst (GWB; Ivan Puttenl (j2001bl ): Ivan Putten &: Levinson 
(120031 )). For stellar mass black holes, this output may be detected by advanced gravitational 
wave detectors for events in the local Universe. As candidate inner engines to long GRBs, 
evidence for the associated spin down of the black hol e has been found in the no rmalized 
light curve of 600 long GRBs in the BATSE catalogue (Ivan Putten fc Guptall20~09l ). 



The high frequency r ange of the planned advanced detectors LIGO- Virgo (IBarish &: Weiss 



1999 



Arcese et alll2004f) . the Large-scale Cryogenic Gr avitational-wave Telescope ( LCGT , 
Kuroda et all fcoioh ) and the Einstein Telescope (ET, Iffild. Chelkowski fc Freisel (boo 



covers the quadrupole emission spectrum of orbital motions around stellar mass black holes, 
thus establishing a window to rigorously probe the inner most workings of GRBs and some 
of the CC-SNe. 

We anticipate an event rate for long GWBs of 0.4-2 per year with in a distance of 100 Mpc 
from the local event rate of long GRBs (IGuetta fc Delia Valid 120071) . I t compares favorably 



with t hat of mergers of binary neutron stars (e.g lO'Shaughnessy et all (120081 ): lAbadie et al 
(120101 )). Since the local event rate of type Ib/c supernovae is ~ 80 per year within a distance 



of 100 Mpc, the branching ratio of Type Ib/c supernovae into long GRBs is therefore rather 
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small, about 0.5 (Ivan Puttenl 120041 ) up to ~ 2.5% flGuetta fc Delia Valid 120071 ) . It suggests 
the existence of many failed GRB-supernovae, notably supernovae with relativisti c ejecta, su- 



pernovae with pronounced aspherical explosions and radio- loud supernovae (e.g. iDella Valle 
(120101 ) and references therein). Therefore, the rate of events of interest to potential bursts 
in gravitational waves appears to be 1-2 orders of magnitude larger than the event rate of 
successful GRB-supernovae. In addition, also type II SNe, whose event rate is 3-4 time s 



1999 



larger than that of type Ib/c (jCappellaro. Evans fc Turatto 
may explode and expand asymmetrically ( iHoflich. Wheeler fc Wang 



Vlannucci et al. 2005). 



1999 



Ishikawa et al. 



19921 ). suggesting a significant additional potential for gravitational- waves burst production. 



A blind rather than a triggered search for bursts events in the local Universe seems to 
be appropriate by taking advantage of the all-sky monitoring capability of the gravitational- 
wave detectors in view of a beaming factor of gamma - ray bursts of fh < 10 ( 9 > 2 5 
deg) up to a few hundred ( 9 ~ 4 deg, iFrail et al.l (120011 ); Ivan Putten fc Regimbaul (120031 ); 
Guetta fc Delia Valid (120071 )). A blind search is also expected to be competitive with current 
X/optical surveys for detecting the shock break-out associated with an emerging CC-SNe, as 
it lasts only a few dozens or minutes up to a few hours, and it naturally includes the possibil- 
ity of long GRBs coming from merger events with no supernova, which may be exemplified 
by the long event GRB060614 of duraton 102 s discovered by Swift. 

While the energy output in long gravitational wave bursts (GWBs) produced by rapidly 
rotating black holes should be large, searching for these bursts by matched filtering is chal- 
lenging in view of anticipated phase-incoherence due to turbulent magnetohydrodynamical 
motions in the inner disk or torus. 

Here, we focus on the detection of a trajectory in the time frequency domain produced 
by long GWBs, satisfying phase-coherence on short up to intermediate timescales. This 
objective goes further than the detection of a burst signal, with the aim to extract reason- 
ably accurate information on the burst evolution. Inevitably, the sensitivity distance for 
extracting trajectories is considerably more conservative than the sensitivity distance for a 
detection per se. 

We shall discuss a new matched filtering detection algorithm to detect trajectories in 
the time frequency domain for long GWBs with slowly varying frequencies lasting tens 
of seconds with intermittent phase coherence. For a burst lasting 50 s, for example, the 
algorithm searches by matched filtering using segmented templates on a time scale of, e.g., 
1 s. This procedure gives a compromise between optimal matched filtering, applicable to 
phase-coherence extending over the entire burst duration as in binary coalescence of two 
black holes, and second order methods by correlation of independent detector signals in the 
time-domain. For our example, the compromise results in a sensitivity distance below that 



-4- 



of optimal matched filtering by a factor of about v^O ~ 7, and an improvement by a factor 
of VlOOO ~30 over second order methods for signals around 1 kHz. 

In Section 2, we discuss the astronomical origin of long GRBs from possible both CC-SNe 
and mergers. In Section 3, we introduce a model and template for long GWBs from rapidly 
rotating Kerr black holes. In Section 4, we describe the proposed time sliced matched filtering 
search algorithm and the evaluation of the sensitivity distance for a reasonably accurate 
extraction of trajectories in the time frequency domain. Our findings are summarized in 
Section 5. 



2. Long GRBs from CC-SNe and mergers 

As a universal inner engine, rapidly rotating Kerr black holes can explain long GRBs 
from both CC-SNe and some of the mergers. They enable long GRBs with supernovae exclu- 
sively in star forming regions with stellar wind host environments and without supernovae 
such as GRB060614, both in and away from star forming regions including the halo such as 
GRB 070125. 

Some recent Swift and Fermi-LAT detections of long GRBs show events with and with- 
out X-ray afterglows such as GRB050911, ISM-like constant or wind-like r~ 2 density profiles 
in the local host environments, and a diversity in distances from a local host galaxy, as 
summarized in Table I. This phenomenology is not accounted for by CC-SNe alone, and sug- 
gests that some of the long GRBs are associated with mergers, otherwise sharing a common 
long-lived inner engine. 

Radiative processes around Kerr black holes are driven by frame-dragging, expressed in 
terms of a non-zero angular velocity u of zero-angular momentum observers. It gives rise 
to high energy emissions along the sp in axis by the induced pote ntial energy E = QJ p for 
particles with angular momentum J p (Ivan Putten &: Gupta 120091 ) and, contemporaneously, 



to a spi n connection to su rrounding matter by equivalence in topology to pulsar magneto- 



spheres (Ivan Putten! Il999l ). Thus, observations on high-energy emissions (in gamma-rays, 
afterglow emissions, possibly ultra-high energy cosmic rays) from leptonic jets as well as on 
low-energy emissions (in gravitational waves, MeV-neutrinos and magnetic winds, powering 
a supernova or a radio-burst) emanating from matter is required for full calorimetry. 

According to the Kerr metric, rotating black holes are energetically very similar to 
spinning tops in the sense of having a ratio 

Erot - - cos- 2 (A/4) = 0.5 - 0.5858 (1) 
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TABLE I. Proposed core-collapse and merger progenitors to a Swift sample of long GRBs. 



GRB 


Redshift 


Duration [s 


host 


constraint" 


type 


050820A 1 ' 2 


1.71 


13 ±2 


UVOT < 1 arcsec 


ISM-like 20 


merger 


050904 3 ' 4 ' 5 


6.29 


225 ±10 


unseen low SFR 


dense molecular cloud 21 


CC-SN 


050911 


0.165 


16 


cluster EDCC 493 


no X-ray afterglow 6 


23 24 

merger ' 


060418 7,8,9,io 


1.490 


(52 ± 1) 


ISM spectrum 


7-ray efficiency 20 


merger 


060505 11 


0.09 


4 


spiral, HII 


no SN 13 


merger 


060614 11 ' 12 ' 13 ' 14 


0.13 


102 


faint SFR 


no SN 12 ' 13 


merger ' 


070125 15 ' 16 


1.55 


> 200 


halo 


ISM-like 22 


merger 


080319B 1 ™ 9 


0.937 


50 


faint dwarf galaxy 


wind 20 


CC-SN 



a ISM-like refers to a cons tant host density, wind re 



with a massi ve proge n itor ( Sari 



(2005). 5 Berger et al 



(2007aJ); 



( 2007). lCummungs et al 



Berger et al 



ersto a r density profile associated 



2000) 



19981: Chevalier fc Li 
Band et all (l2005ft . ffivnbo et aD toO^) : ■Iciisiui > « t al 



(|200J), l 1Covirio et all (120081): 



Cenko et al.l (12 008) 



1 .5 m telesco p e. iFynbo et al.l (12 006 



Bloom et al. 



16 



Chandra et al. 



(2008), 1£ 

( 2008 ): 2 ^Paee et al . 

(l2008h . 2 6 iCaito et al. I tom . 



14 



Gehrels et al. 



tea 



17 



(I2009h : 2C ICenko et all (12010ah. 2] lFrail et a 



(2007), 12 observed with a 8.2 m tele scop e, Delia Valle et al. ( 2006h. 13 observed with a 



ii 



( 2006) . 4 Sak amoto et al. 



d2007bl) : ^Fa lcone et al.l (Eooii .iProchaska et al.l 



Jakobsson fc Fynbo 



(2006]), 1 



Racusin et a. 



(120061 ). 24 lvan Putten fc Gupta! (120091 ): 



Gal- Yam et al.l (120 06): 



fboosl) ^Tanvir et al. 



fl2006l). 2 1Chandra et al. 



van Putten 
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close to the Newtonian value of | for all spin rates Qh associated with the angular momentum 
J and rapidity A, where sin A = a/M, a = J/M. In the absence of a small parameter, E rot 
reaches 29% of the total mass-energy of an extreme Kerr black hole, or about 1.6 — 6 x 10 54 
erg for a black hole in the mass range 4 — 14M Q . 

Existing calorimetry on GRB-afterglows points to true energies in gamma-rays E 1 and 
kinetic energies Eke broadly distributed around 1 x 10 51 . These estimates apply to the 
sub-sample of GRBs for which achromatic breaks in their light curves have been deter- 
mined, as the basis for estimating the half-opening angle of the beamed outflow by a 
mod el fit for the intera ction with the host environment comprising about seven parame- 
ters (jCenko et al.ll2010af ). It has been well-recognized that there is considerable uncertainty 
in this estimation procedure, also in assuming a uniform luminosity profile across the out- 
flow. Quite generally, the outflow may comprise a baryon-poor ultra- relat ivistic inner jet 



and a possibly baryon-ric h and mildly rela ti vistic collimating wind fe.g. iPederson et al. 



rll998h: [Frail et all (boool): kumar fe Piranl (hood ) ; [Ramirez- Ruiz. Celotti fc Reesl (l2002h: 
van Putten fc Levinsonl (120031 ) ; iBerger et al.l ( 120031 ) i lHuang et al.l ( 120041 ) ; |Peng. Konigl fc Granot 
( 120051 )). while the inner jet may be structured in different wa ys, e.g., with the highest lu 



minosity reached at the boundary with the collimating wind ( van Putten Sz Gupta 2009h . 



Despite these uncertainties and model assumptions, exceptional events (jCenko et al.ll2010af ) 



point to £ 7 ~ 2.25 x 10 52 erg (GRB060418) and E KE ~ 3.56 x 10 52 erg (GRB050820A), well 
above their typical values around 1 x 10 51 erg. And yet, these energies represent a minor 
fraction of about 1% of E rot or less following (CQ). 

The common properties in durations, spectra and true energies in the prompt gamma- 
ray emissions in long GRBs point to a universal inner engine as the outcome of the apparent 
diversity in astronomical origin and host environment (Table I), whose lifetime is intrinsic 
to the physical state of the energy reservoir itself. By the above, identifying it with the 
rapid spin of the central black hole leaves about 99% of the spin energy unaccounted for 
and a similar fraction "unseen" in emissions in gravitational waves and MeV-neutrinos, the 
remainder of which may be dissipated in the event horizon. 



3. Long duration gravitational-wave bursts from Kerr black holes 



The lowest order, axisymmetric component o f magnetic fields in an accretion disk can 
form a torus mag n etosphere around a bla c k hole (|van Puttenll2001bl). In its low est energy 



state (IWaldl 1 19741 iRuffini & Wilson! 119751 : iDokuchaevI 119871 Ivan Putted l2001bj ). a rapidly 
spinning black hole can hereby establish a spin connection to its surrounding disk or torus 
for catalytic conversion of a major fraction of its spin energy into a variety of radiation 
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channels. 



Quite generally, the torus magnetosphere will be highly time variable. When the central 
black hole spins rapidly, the interaction is governed the variance of the net poloidal magentic 
flux of the magnetic field to mediate energy and angular momentum transfer. The black hole 
hereby sweeps up the inner face of the torus to excite non-axisymmetric instabilities, giving 
rise to gravitational radiation, MeV-neutrinos and magnetic winds. This state of suspended 
accretion lasts for the lifetime of rapid spin of the black hole, and sets the duration of a 
possibly accompanying GRB. During viscous spin down, a major fraction if not most of the 
spin energy is dissipated unseen in the event horizon, creating a astronomical amounts of 
Bekenstein entropy. 

Based on the Kerr metric, the model template produced by matte r swept up by a ce ntral 
black hole shows negative chirp with late time asymptotic frequency ( Ivan Puttenl 120081 ) 



GW 



590 Hz ( 



10M P 







M 



(2) 



where M denotes the initial mass of an initially rapidly spinning black hole. The negative 
chirp represents the expansion of the inner disk or torus during spin down of the black 
hole on a viscous time scale. Thus, the merger of two neutron stars features a late-time 
frequency of 1.5-2 kHz defined by the sum of the mass of the two pr ogenitor neutron stars 
(2 x 1.5 - 2M Q , where the high mass es refer to PSR J0751+1807 JNice et all l2005h and 
PSR J1614-2230 ( faemorest et aJboiol )). Lower asymptotic frequencies down to about 500 
Hz are produced by high mass black holes formed in CC-SNe, or mergers of neutron stars 
with a high mass black hole companion. Mergers and CC-SNe hereby represent distinct 
astronomical progenitors of long GWBs (and of long GWBs) with a diversity in asymptotic 
frequencies corresponding to a range in the black hole mass. 

The full template follows by integration of the equations of conservation of energy and 
angular momentum i n the process of b lack hole spin down mostly against the inner face of 
a surrounding torus (Ivan Puttenl 120081 ) 



m = -k(q h -n T )n T , J = -k(q h -n T ), 



(3) 



where k represents the spin connection due to the (variance) in poloidal magnetic field for 
a black hole in its lowest energy state and the angular velocit y fir of the toru s is tightly 
correlated to that at the ISCO. As non-axisymmetries develop (Ivan Puttenl 120081 ) . M and J 
will be carried off mostly in gravitational waves, until the fixed point Qjj — &>t is reached 
at low a/M. The resulting gravitational- wave templates obtain by integration of ([3]) as a 
system of two coupled ordinary differential equations in response to an inital mass and spin 
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rate. In Ivan Puttenl (120081 ). we calculated the model templates for different initial values 
of a/M. The resulting strain amplitude and frequency as seen at Earth for a source at a 
fiducial distance D — 10 Mpc produced by a black hole with initially maximal spin (a = M) 
is shown in Fig. 1. Here, we show the orientation averaged strain amplitude, which is a 
factor of 1 / \/5 smaller than the am plitude at optimal orientation of the source along the line 
of sight ([Flanagan fc Hughes! 1 19981 1 . 



In general, the total energy output and the frequency scale with M and M _1 , respec- 
tively. The strain amplitude scales with casa function of the mass of the torus relative to 
M (typically about 0.1-1%), while the total duration of the burst scales with M and K~ l . 



4. Time sliced matched filtering 

The application of matched filtering depends crucially on phase coherence in the true 
signal, in correlating it to a model template. For a magnetohydrodynamical system powered 
by a Kerr black hole, turbulence in the inner disk or torus inevitably creates phase incoher- 
ence on long time scales, inhibiting the application of matched filtering by correlating to a 
complete wave form template with the detector output. To circumvent this limitation, we 
first slice a model template into iV = T 90 /t segments on intermediate time scales, r, for 
which phase- coherence may be sustained. Matched filtering is now applied using each slice, 
by correlating each template slices i with the detector output with arbitray offset S in time. 

We here re port on a test of our alg orithm on strain noise amplitude data of the TAMA 



300m detector ( iTakahashi et al.l l2004j ). TAMA 300m was the first laser interferometric 
gravitational-wave detector with long duration continouus operation at a reasoble sensi- 
tivity of better than 10~ 20 Hz _1//2 in strain amplitude noise. Here, we use data from run DT8 
(2/2003-4/2003) and DT9 (11/2003-1/2004) (partially coincident with the LIGO S3 run). 
The data are sampled at 20 kHz and organized in frames of 2 20 samples (about 52 seconds). 
During DT8-9, we note that some 51 CC-SNe have been observed in the local Universe with 
a mean redshift of 0.024 (D = 100 Mpc). Fig. 2 shows the detector strain noise amplitude 
for DT9. Relevant to our exploration of the detector sensitivity to the long GWBs is the 
high frequency range between a few hundred and a few thousand Hz. Given this focus and 
a noticeable increase in the detector strain noise amplitude below 550 Hz, we apply a band 
filter between 650 Hz and 4000 Hz to the detector strain noise amplitude, where the noise 
is Gaussian to high precision as shown in Fig. [2j To excellent approximation, therefore, the 
detector output h n {t) is representing by the sum of a signal s(t) at the detector and white 
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time [s] time [s] time [s] 



Fig. 1. — (Left.) Shown is the evolution of the orient at ion- averaged strain amplitude of the 
model template of a long GWB produced by a Kerr black hole of mass 10M Q at a fiducial 
distance of 10 Mpc (left top) and the negative chirp (righttop) with a decay in frequency 
(right top) to (J2J) associated expansion of the ISCO up to the fixed point Qh = of ©• 
For time sliced matched filtering, the template is partioned into segments of intermediate 
duration, here sub-windows of 1 s duration of which a few are highlighted by the solid bars. 
Each segment has an effectively band-limited frequency spectrum, shown here for some of 
the 52 segments (middle). The maximum luminosity is reached when a/M ~ 0.8. For 
a/M > 0.8, the luminosity is less than maximal as — ^isco remains close to f2#. For 
a/M < 0.8, the luminosity decays, as black hole continues to relax to a slowly spinning, nearly 
Schwarzschild black hole corresponding to the fixed point Qh — ^r- (Right.) Snapshots of 
the model templates for h n (t) in segments 1, 10, 20 and 40 out of a total of 52. 
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additive noise w(t), 

h(t) = s(t)+w(t), 8 {t)=(^^)s{t), (4) 

where D denotes the source distance and S(t) denotes the template S(t) (0 < t < T 90 ) of 
the signal for a source at a reference distance of 10 Mpc. 

For illustrative purposes, we consider a time sliced matched filtering (TSMF) with in- 
tervals r = 1 s and a burst duration T90 of 50 s. Slicing the template S(t) (0 < t < T90) into 
segments of duration r gives N s = T 90 /r = 50 templates Si(t) for each slice i, 

S i (t) = S(t* + t), t* = ir, (i = l,2,...,N„ 0<t<r). (5) 

Here, it is understood that the time t is discrete, i.e., t = tj, tj = (J/Nf) s, where Nf = 20000 
denotes the 20 kHz sampling rate of the TAMA 300 detector. TSMF over segments of all 
data points (N = N s x Nf ~ 10 6 ) in the TAMA 300 data frames starts with computing the 
Pearson coefficients 

f " (S) = IKlTlNl (6) 

as a function of the unknown offset 5 = z/20000 (i — 0,1, ... , N) over a complete frame of N 
samples obtained at the sampling rate of 20 kHz, representing the uncertainty in continuation 
of phase between the time slices, where 

{St ■ h) s = E 0<tj<T Si(t* + t,)h(5 + tj) (7) 

denotes the discrete inner product of the Si(t* + tj) and the h(5 + tj) with L 2 norms 

115,1 1 = [£o<t,<TS*(*:+ti)] 1/2 , \\h e \\ = [X 0<tj<T h 2 (5 + tj)] 1/2 . (8) 
Thus, pi (5) represents the cosine between the Si and h in (jSJ). 

To define a signal to noise ratio SNRi for each slice i, we proceed with the normalized 
Pearson coefficients 

pi{o) = , (9) 

where cij denotes the standard deviation of Pi(S) as a function of 5, and define 

SNRi = max Pi (5). (10) 

<5 




-5 5 -5 5 

normalized h normalized h 

Fig. 2. — {Top.) Shown is the TAMA 300 strain noise amplitude in the one-sided frequency 
domain in the DT9 run of 2003-4 produced by 1 frame of 52 s duration. (Middle.) The 
distribution is non-Gaussian largely due to low frequency noise below a few hundred Hz. 
(Bottom.) The fluctuations of the observed probability distribution of the strain noise am- 
plitude over 1 second, normalized to a = 1 and plotted over 250 bins of normalized strain 
amplitude noise h n , about to the Gaussian distribution with o = 1 can be compared with 
the expected fluctuations on the basis of the standard error in the mean in each bin (SEM, 
smooth 1-cr curves). A band pass filter across 650 Hz - 4000 Hz recovers essentially Gaussian 
noise. 
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In order to extract a trajectory in the time frequency domain, we proceed with corre- 
lations pi{8) > 4, and select the largest correlation in each time slice. This filter extracts at 
most 52 data points from each frame of 52 seconds when r = 1. 

Fig. 3 shows the statistical results over 10 frames for spin down of a black hole of 
M = 10M© at distances of D = 0.050 Mpc and D = 0.070 Mpc. Here, extraction of a 
time-frequency trajectory is performed by polynomial interpolation of the aforementioned 
< 52 filtered points with rejection of points with large scatter. Fig. 4 indicates a sensitivity 
distance of about D ~ 0.07 — 0.1 Mpc in the black hole mass range M = 10 — 12M . 



5. Conclusions and outlook 



The diversity in multi-wavelength phenomenology on long GRBs strongly suggests a 
common inner engine that is intrinsically long-lived representing the outcome of various 
astronomical scenarios. We here identify the inner engine with rapidly rotating Kerr black 
holes, whose lifetime is set by the secular time scale of spin in a process of spin down 
against surrounding high density matter. At present, the most quantitative observational 
evidence for black hole spin down is fo und by matched filtering an alysis of 600 light curves 
of long GRBs in the BATSE catalogue (Ivan Putten &: Guptall2009l ). This mechanism points 
to major contemporaneous emissions in gravitational waves and MeV-neutron emissions, 
that is likely to be dominant over the presently observed electromagnetic radiation in GRB- 
afterglow emissions and kinetic energies in GRB-SNe. True calimetry on the inner engine 
requires observations in these non-electromagnetic windows. 



Most GRBs remain unobserved, due to beaming factors of "a few" to about 500 (IFrail et al 



200ll ; lvan Putten &: Regimbaull2003l ). For CC-SNe, there is a considerable uncertainty of up 
to several days for the time of onset based on extrapolating backwards in time the super- 
nova optical l ight curve. Not all long GRBs appear to be associated with CC-SNe, i.e., 



GRB 060614 flDella Valle et al.ll2006l : iFvnbo et al. 



2006; Gehrels et al. 2006; Gal- Yam et al. 



20061 ) appeared without a bright supernova. If GRB 060614 was not an anomalously fain t 
("failed") supernova, it may have been a merger event (Ivan Puttenll2008l ; ICaito et al. II2009I ). 
Similar considerations apply to the halo event GRB070125, which ap peared with no opti- 
cally identified host galaxy ( Cenko et al.l 120081 ; IChandra et al.l 120081 ). A significant frac- 
tion of binar ies exist in globul a r clus ters as indicated by their population of luminous X- 
ray sources (jVerbunt &: Lewinl (120061 ) and references therein) . As the number of globular 



clusters correlates with th e luminosity of elliptical galaxies ([Harris fc van den Bergi 11981 



Burkert fc Tremaind |2010| ). the latter in particular may be preferred sites for GRB070125 
type events. Therefore, blind, untriggered searches appear to be preferred in searches for 
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Fig. 3. — {Top.) Time frequency trajectories in the time-frequency domain extracted by time 
sliced matched filtering (TSMF) for spin down of a M = 10M Q black hole to a remnant of 
8AM & . The results shown in the top windows represent the outcome of injections into 10 
TAMA 300 frames of 52 second durations for two source distances. (Middle.) The extraction 
is based on the largest SNR> 4 in Pi(5) (on < 5 < 23 s, containing most of the signal) and 
at most one data point in each of the r = 1 s bins i = 1, ... ,51 (vertical lines), followed by a 
polynomial interpolation (dotted curve) of points with reduced scatter (circles). (Bottom.) 
The signal to noise ratios of the data points used in the polynomial fit. 
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0.05 0.055 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 0.1 



distance [Mpc] 

Fig. 4. — Shown is the resulting deviation of the extracted time-frequency trajectory from 
the model trajectory of the injected signal as a function of distance D = 0.05 — O.lMpc 
and black hole masses M = 8 — 12M . The curves represent the standard deviation of the 
deviations in the extracted trajectories over a common time interval 12 — 16 s. The results 
show sensitivity distances from D = 0.05 Mpc for M = 8M & to about D = 0.1 Mpc for 
M = 12M and the average trend line (thick line). 
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TABLE II. Estimated 1 — o uncertainty in the extracted time frequency trajectories by 
TSMF as a function of distance applied to long bursts in GWs produced by black hole spin 
down against high density matter, expected to form in some of the CC-SNe and mergers of 
neutron stars with a rapidly rotating companion black hole. We do not include results on 
neutron star-neutron star mergers, in view of their relatively high frequencies of 1.5-2 kHz 
away from the region of maximal sensitivity of the existing gravitational wave detectors. 



Mass (M ) 


£>[TAMA] a (Mpc) 


D[Adv] 6 (Mpc) 


R C D 


a 


max(SNRi) d 


sum(SNRi) e 


8 


0.05 


25 


0.1 


6% 


6.4 


74 


8 


0.07 


35 


0.3 


22% 


5.0 


54 


8 


0.10 


50 


1.2 


27% 


4.5 


46 


10 


0.05 


25 


0.1 


4% 


8.2 


96 


10 


0.07 


35 


0.3 


11% 


6.2 


68 


10 


0.10 


50 


1.2 


22% 


4.9 


43 


12 


0.05 


25 


0.1 


<1% 


10.5 


130 


12 


0.07 


35 


0.3 


4% 


7.8 


87 


12 


0.10 


50 


1.2 


9% 


6.0 


61 



a With h n ~ 10~ 21 Hz~5 at 1 kHz during DT8 (2/2003-4/2003). 
b With h n = 2x 10~ 24 Hz~ 5 at 1 kHz. 

c Estimated event rate within distance _D[Adv] assuming 10 times more relativistic CC-SNe 

than successful GRB-SNe with otherwise the similar inner engines and the observed event 

rate of 1 long GRB per year within D = 100 Mpc. 

d Based on Pi(S), < 5 < 23 s, r = 1 s and averages over 10 frames. 

e Based on SNRj > 4 and averages over 10 frames. 
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local events, exploiting the all-sky survey capability of gravitational wave detectors. 

To explore the sensitivity distance of advanced gravitational-wave detectors to the an- 
ticipated long duration negative chirps from events in the local Universe, we introduce a time 
sliced matched filtering algorithm and apply it to the strain noise amplitude of the TAMA 
300 detector with signal injection. A time sliced approach can circumvent the limitations 
posed by phase- incoherence on the time scale of the duration of tbe burst, which generally 
inhibits optimal matched filtering using complete wave form templates. By injecting our 
model template into the strain noise data of the TAMA 300 detector during a run with 
h n {t) ~ 1 x 10~ 21 Hz~5 ; we compute correlations between the i— th time slice (of duration 
r) and one frame (about 52 s) of detector output. 

For a typical black hole mass of 1OM , the results indicate the importance of the middle 
and late time behavior of the bursts (Fig. 3) below 2000 Hz, but less so the initial spin 
down phase associated with higher frequencies when starting from a maximal rotation. This 
reduction results from the relatively large strain amplitude noise of the detector in the shot- 
noise region for frequencies of a few kHz. Effectively, the search is focused on the output 
post-maximum (a/M < 0.8), which obviates the need to consider templates over a wide 
range of initial spin. A complete search will include a scan over two parameters: a range 
of intermediate time scales, e.g., r = 0.1 — 1, and scaling of frequencies of the templates to 
account for a diversity in black hole masses, e.g., between 5-15 M & . The range r = 0.1 — 1 
s represents tens to hundreds of orbital periods of the torus. This time scale in coherent 
evol ution of the t o rus m ight correspond to that of sub-bursts in the light curves of long GRBs 
(e.g. Ivan Puttenl (119991 )). and quasi-steady evolution of the torus on a time scale of at least 



tens of orbital periods follows from an upper bound of about 10% on the electromag netic 



field energy it can support relative to its kinetic energy (Ivan Putten fc Levinsonl 120031 ). 



Our estimated TAMA 300 sensitivity distance for extracting time frequency trajectories 
is summarized in Table II. The results show D ~ 0.070 — 0.10 Mpc for black holes in the mass 
range M = 10 — 12M Q , which compares favorably with the sensitivity distance for neutron 
star-neutron star coalescence. Extrapolation points to a sensitivity distance D ~ 35—50 Mpc 
for a strain noise amplitude of 2 x 10~ 24 Hz~5 in the planned advanced detectors LIGO- Virgo 
and the LCGT. This sensitivity distance serves as a conservative estimate for the sensitivity 
distance for a detection. In particular, the sums of the 15 SNRi's shown in Fig. 4 are 95 and 
77 for D = 0.05 Mpc and, respectively, D = 0.07 Mpc. These sums point to a sensitivity 
distance for a detection (with no particular information of behavior in the time frequency 
domain) on the order of a few tenths of Mpc, corresponding to sensitivity distance of a well 
over 100 Mpc for the advanced detectors. 



Based on the observed event rate of long GRBs of about 1 per year within a distance 
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of 100 Mpc, the observable event rate suitable for extracting time frequency trajectories will 
depend on the abundance of the parent population of aspherical, relativistic and radio-loud 
Type Ib/c that may be powered by irrad iation of the stellar envelope from a long-lived 
black hole inner engine (Ivan Puttenl 120031 ). It may reach one per few years if their event 
rate is about one order of magnitude larger than the rate of successful GRB-supernovae. 
The observable event rate could be larger if a fraction of the supernovae of Type II is 
similarly powered by long-lived black hole inner engines. A much larger sensitivity distance 
is anticipated for the planned 10 km ET in Europe. Conceivably, all sky radio surveys, e.g., 

will furthe r provide us with a probe for long 



the LOw Frequency ARray (ILOFAR 



duration radio bursts from mergers ( 



2010) 



van Puttenl 120091 ) and, combined with gravitational 



wave surveys, provide a direct measurement of the relative event rate long GRBs from 
mergers to long GRBs from CC-SNe. 



The diversity in the origin of long GRBs in both CC-SNe and mergers (Ivan Puttenll2009l ) 
and the comparible sensitivity distances of their potential emissions in gravitational waves 
and those from binary coalescence suggests the need for extended searches over the complete 
frequency range of both. Apart from scaling by black hole mass and a diversity in initial 
spin, the proposed long GWBs are universal, and their progenitors are revealed only by the 
absence or presence of a precursor signal in gravitational waves in case of, respectively, a 
CC-SNe or merger event. 
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